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We discuss prospects for direct measurement of stochastic gravitational wave background around 
0.1-lHz with future space missions. It is assumed to use correlation analysis technique with the 
optimal TDI variables for two sets of LISA-type interferometers. The signal to noise for detection 
of the background and the estimation errors for its basic parameters (amplitude, spectral index) are 
evaluated for proposed missions. 



I. INTRODUCTION 



Gravitational wave background generated in the early universe is one of the most fascinating targets in observational 
cosmology ,1, ^]. Among others, inflationary theory has a realistic mechanism to generate the background, and indeed 
confirmation of the background is regarded as another strong support for presence of inflationary phase in the early 
universe. While we might indirectly detect the inflationary generated background by B-mode polarization analysis 
of CMB ill, direct detection of the background with gravitational wave detectors is an indispensable approach to 
study the inflation in more detail. This is because the amplitude of the background is mainly determined by the 
value of the inflation potential when the gravitational waves cross the Hubble horizon in inflationary epoch. If we can 
measure the amplitudes at two widely separated frequencies {e.g. ^ 10^ ^^Hz and ^ O.lHz), the global structure of 
the potential might be constrained 0,ISI3 (see also '»]). The slope of the spectrum at a given band will also provide 
us information of the derivative of the potential. Therefore, it is quite meaningful to understand how well we can 
measure the basic parameters that characterize the spectrum. 

Standard slow-roll inflation predicts that the spectrum ^Gw{f ) is nearly flat at frequency regime relevant for the 
direct detection ([S, '6], see 9, 10, 11, 12] for recent studies). This means that its strain amplitude is expected to be 
higher at lower frequencies. However, astrophysical foreground would be a fundamental obstacle to directly detect 
weak inflationary background below ~ O.lHz. For this reason the band around ~ IHz is considered to be suitable 
for the direct detection, and projects such as, the Big Bang Observer (BBO; US) or DECI-hertz interferometer 
Gravitational wave Observatory (DECIGO; Japan) have been proposed (see also 0,^^). For these projects, 
correlation analysis is a powerful method to observe weak background 0, 0, 0, . In this paper prospects of this 
method are studied quantitatively, from detectability of the background to parameter estimation errors. 

This paper is organized as follows; In section II we study basic aspects of the optimal data streams for the Time- 
Delay-Interferometry (TDI) method, and designed sensitivities of BBO or DECIGO are briefly mentioned. In section 

III a formal discussion for correlation analysis is presented. We evaluate the expected signal to noise ratio for detecting 
the background and derive expressions for parameter estimation errors based on the Fisher matrix approach. In section 

IV numerical results for BBO project are given with using formulas in section III. Section V is a brief summary of 
this paper. 



II. TDI VARIABLES AND THEIR RESPONSES TO GRAVITATIONAL WAVES 

First we summarize standard notations to discuss stochastic gravitational wave background po|. The plane wave 
expansion for gravitational waves is given by 

K,{t,x)= r df f dnhA{L^)e'"^^'-^-''^e^,{n), (2.1) 

p=+ X s'^ 

where / is the frequency of each mode, fl is the unit vector for its propagation, is a unit sphere for the angular 
integral dft, and e^^ (P = +, x) is the basis for the polarization tensor. We assume that the stochastic background 
is isotropic, unpolarized and static, and express the spectrum of its amplitude /i^i (/,$!) in terms of the logarithmic 
energy density of the gravitational waves ^cwif) = ^/PcdpGw{f )/dlnf {p^. the critical density) as follows; 

{h*p{f,mp'{f\^')) - ^s\n-n')Spp,6{f-n\f\-^nGw{\f\), (2.2) 
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where Hq is the Hubble parameter and we fix it at Hq =70km/sec/Mpc. The symbol (...) represents to take an 
ensemble average of stochastic quantities. 

Next we discuss responses of interferometers to incident gravitational waves. We concentrate on a LISA- like detector. 
Each unit is formed by three spacecrafts at the vertices of a nearly regular triangle, as shown in the solid lines in 
figure 1 where we also define the labels of its vertices (1,2,3) and arms {Li, L2, L^). We denote the six (one-way) 
relative frequency fluctuations of the laser light as yij{t) {i, j — l,2,3;i ^ j). The quantity yij{t) corresponds to the 
signal measured at the spacecraft j, transmitted from the spacecraft i, 7^ j) along the arm Li. For example, the 
variable 2/13 responds to a single gravitational wave mode with parameters (/, ft and e^^) as 



1 ni2 • ■ ni2 
21-(JI-ni2)2 



(1 + nni2)/iA(/, mU{t. 1) - U{t - T, 2)), 



(2.3) 



where ni2 is the unit vector from vertex 1 to vertex 2, and the function U{t,i) contains information of the phase of 
the wave at time t and position Xi of vertex i. It is given as 



(2.4) 



The time delay interferometry (TDI) is an important technique for LISA-type detectors to overcome the laser frequency 
fluctuations j23|. We follow Ref . |23l| to summarize the relevant data streams for signal analysis (see also |24L l25L l2fl| V 
We first define a TDI variable, a as follows; 



a = 2/21 - 2/31 + 2/13,2 - 2/12,3 + 2/32,12 - 2/23,13, 



(2.5) 



where we used the notations like 2/13.2 (^) = Viait — L2) and 2/32,12 (i) = 2/32 (^ — Li — ^2)- In the same manner we define 
two other TDI variables (3 and 7 that are given by cyclic permutations of the subscripts of the variable a. These a, 
/3 and 7 are TDI variables, but their noises are correlated. Thus we define the following new variables A, E and T as 



A 



E = _(a-2/3 + 7), 
T = -L(« + ^ + ^). 



(2.6) 
(2.7) 
(2.8) 



As we can easily confirm with using the symmetry of the original data streams a, P and 7, the noises of the variables 
A, E and T do not have correlation. In other words, they are orthogonal. We regard them as the fundamental data 
sets for correlation analysis. Hereafter, we do not discuss differences of the arm lengths Lj or their time variations, 
and simply put Li — L2 — L3 = L ~ const. Even if the second generation TDI variables [23, 123 are used, our basic 
results are not changed 24] . As our primary interest is observation of the monopole mode of the background, the 
motion of the triangle will not be included. 

At low frequency regime with fL <C 1, the response of E mode is approximately given as |27| 



E'- 



i(nx ■ h ■ nx - ny ■ h ■ ny) 



(2.9) 



where directions of two unit vectors nx and ny (nx -L ny) are given in flgure 1. The large parenthesis in ea. (|2.9|) 
is the response of a simple i-shaped detector, as often used in the literature of gravitational waves. In the same 
manner, the A mode asymptotically becomes a simple response obtained by rotating nx — ny system by 45° on the 
detector plane :30j. The response of T mode at low frequency regime is T = 0{f^) in contrast to A = E ~ 0{p) as 
in ea. H2.9|l . This allows us to use the T mode to monitor the detector noise in principle [sills^. 

As for the sources of noises, we take into account the proof mass and optical path noises with parameters Wp and 
Wn as 



gproof —mass 
goptical—path 



2.5 X 10-''^ii;2(//lHz)-2Hz-\ 
1.8 X 10"^^w^(//1Hz)2Hz~^ 



(2.10) 
(2.11) 



The proof mass noise is dominant at low frequency regime. Basic parameters of LISA are {L,Wp,Wo) = (5.0 x 
lO^km, 1.0, 1.0) [2^, 27, 33]. For BBO project three possible configurations are discussed 0, namely BBO-lite: (2.0 x 
lO^km, 0.1, 4.0 X 10-*^), BBO-standard: (5.0 x lO^'km, 0.01, 1.4 x lO^^) and BBO-grand (2.0 x loW, 0.001, 2.6 x 10"^) 



FIG. 1: Three vertices (1,2,3) and three arms [Li, L2, L^) of the first unit (solid fines). Configuration of second unit (short- 
dashed fines) is obtained by 180° rotation of the first one around its center. Labels for vertices and arms of the second ones 
are transported with this rotation. 



0. These parameters for BBO configurations were presented in 2003, and we use them for a reference. We should 
keep in mind that they do not always represent the most up to date designed sensitivities. In the original proposal 
the second configuration; "BBO-standard" is simply named as "BBO" . In this paper we use the former for the name 
of a specific configuration and the latter for the name of the project itself. The parameter wq is assumed to scale as 
X^^^{epY^/'^LD-^ with the l aser power p at wavelength A, optical efficiency e, the armlength L and the mirror 
diameter D M. The noises (I2.10[l and (|2.11|) are defined for the one-way signal yij. For the A, E and T modes the 
noise spectra are given by 23, .24., 27] 

SAif) - SeU) = 16sin2(//2/,) [3 + 2cos(///.) + cos(2///,)] 5^™°/-'"--(/) 

+8sin2(//2/,) [2 + cos(///,)] 57*-'^'-f'^*''(/), (2.12) 
Srif) - 2 [1 + 2cos(///,)]' [4sin2(//2/,)5f °°^-'"''^^(/) + 5°P*-'^'-f"*''(/)] (2.13) 

with = (27rL)~i. We have /* = 0.95Hz for BBO-standard and /, = 2.4Hz for BBO-grand and BBO-lite. 

As the responses Hj (/ = A, E, T) to the coefficient /i(/, Vl) of each gravitational wave mode are linear, we can 
express them in a form 

Hi = hp{f,n)R{IJ,n,P). (2.14) 

After taking the ave rag e with respect to the direction and polarization P of the incident waves, we obtain the 
effective noise curve [23 0, 

where the symbol {■■■)^ represents the above mentioned average. Due to symmetry of the data streams {A,E,T) we 
have (i?*(/, /, ft, P)R{I, /, rj, P))^ = for / ^ J. This means that the correlation between {A, E), {E, T) and (T, A) 
vanish for isotropic component of the background. Thus the optimal sensitivity to the isotropic gravitational wave 
background with all the three variables becomes 

In figure 2 we show some of the results for the three possible BBO configurations (see also [l3|). The A and E modes 
have the same sensitivity. In the low frequency regime f «C /* the contribution of the T-mode is negligible to the 
optimal sensitivity, and we have hopt{f) — hA{f)/^ [23 ■ The currently designed sensitivity of Fabry-Perot type 
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FIG. 2: The optimal sensitivities for BBO-lite (thin solid curve), BBO-standard (solid curve) and BBO-grand (thick solid 
curve) configurations. The sensitivities for (A,E) modes (long-dashed curve) and T mode (short-dashed curve) are also given 
for BBO-standard configuration. 

DECIGO is similar to that of BBO-standard [s^. While the minimum noise floor of DECIGO extends to a higher 
frequency (~7.5Hz) than BBO-standard (see figure 2), this difference is not important for observing the inflationary 
background and we will have similar results for these two cases (in order of magnitude sense). Therefore we do not 
take up DECIGO separately from BBO-standard. 



III. CORRELATION ANALYSIS WITH STAR-LIKE CONFIGURATIONS 

Correlation analysis is a powerful approach to detect weak stochastic gravitational wave background ITM IT^ l20| . 
If noises of two data streams have no correlation but their responses to gravitational wave background are correlated, 
we can increase the signal to noise ratio for detection of background by a long term observation. In the case of LISA 
the orthogonal data streams A, E and T do not have correlated responses as mentioned in the last section, and we 
cannot perform correlation analysis to measure the monopole mode of the background. The situation is different 
for studying anisotropies of the background, and the correlation between two of the three {A, E, T) modes will be 
useful to extract information of the background generated by Galactic binaries UMM 

For example, we can show 

that correlation between A and E modes has sensitivity to the hexadecapole mode {I = 4) under the low frequency 
approximation. By combining this effect with the time modulation of the signal due to the rotation of the detector, 
we can get information of low-^ moments ( e.g . I = 2,4) of the Galactic gravitational wave background around ImHz 
that is expected to be highly anisotropic |37ll3a . l39l| . 

For BBO it is proposed to use another unit for correlation analysis to detect isotropic mode. In figure 1 we show 
the proposed placement of its two units 0, 0, 113 ■ Two units have identical specification. The position of the second 
unit (dotted lines) is obtained by rotating the first one by 180° degree with respect to the center of the triangle. We 
transport the labels of the first one to the second one (1', 2', 3', L[, L2, L3) with this rotation. For example, the vertex 
2' is at the opposite side of the vertex 2 around the center. The orthogonal TDI variables A' , E' and T' of the second 
unit are defined in the same manner as the first one in the last section. 

To discuss the correlated response of two variables / and J, we define the overlap reduction function 7//(/) as 

[Ullllsl 

jjjif)^5{R}Rj)^. (3.1) 

Due to symmetry of the relevant data streams, we have 5 {R*iRj)a = 5 (RiRj)^ and the functions 7/j take real 
numbers. Furthermore, their nonvanishing combinations 7/j (/ 7^ J) to isotropic mode are only jaa' ~ Iee' and 
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FIG. 3: The normalized overlap reduction function 7(/) for a star-like constellation. These curves have a scaling parameter 
/* = c/{2nL) that is determined by the arm-length. 

"fTT'- Other combinations {e.g. ^ae'i Ita') become zero. As in the case of LISA, we can, in principle, study 
anisotropic components [l > 1) by correlating two outputs, such as, {A,E) or {A,E'). However, it would be difficult 
to observe the cosmological anisotropy of the background around ~ IHz, conside ring its expected magnitude ~ 10~^ 
[Sflj l (see also [40]). The factor 5 in ea. (|3.1(l is the conventional choice 0, I^J, and we have 77 j = 1 for two 
co-aligned detectors at low frequency limit with the simple response function in the large parenthesis of ea. H2.9|l . In 
figure 3 we show these three non-vanishing overlap reduction functions normalized by 9(///*)'*/2 that comes for the 
prefactor in ea. H2.9|l . These functions are determined purely by the geometry of the placement and scale with the 
parameter = 1/2'kL for any star-like constellation given in figure 1. 

Now we study correlation analysis in a more detailed manner with Fourier space representation. Each data stream 
Si if) is made by gravitational wave signal Hi{f) and noise ni{f) as 

si{f)^Hi{f)+ni{f). (3.2) 

The noise spectrum 

{n*Af)ni{n) = \5{f^nSi{f) (3.3) 

is given in eas. (|2.12|l and (|2.13(l for the orthogonal TDI variables / = {A, E, T, A' , E' , T'}. We assume that the noises 
have no correlation (namely {ti*jTIj) — for / ^ J), and the amplitude of the signal {Hj{f)* Hj{f)) is much smaller 
than that of the noise {nj{f)*nj{f)). The latter is the condition where correlation analysis becomes very powerful. 
We divide the positive Fourier space into frequency segments Fi [i = 1,...,N) with their center frequencies {fi} 
and widths {Sfi}. In each segment the width 6fi is much smaller than fi, and the relevant quantities (e.g. ^cwif): 
luif)) ^re almost constant. But the width is much larger than the frequency resolution A/ = T^^^ (Tabs'- observation 
period) so that each segment contains Fourier modes as many as 6fi/Af ^ 1. 

For correlation analysis we compress the observational data s/(/) by summing up the products S/(/)s,/(/) (I 7^ J) 
in each segment Fi as 

= E (3-4) 

where we omitted the apparent subscript {IJ} for the compressed data {fii} for notational simplicity. As the noises are 
assumed to be uncorrelated, the statistical mean (/i^) is arisen by gravitational wave signal. After some calculations 
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with using eas. l|2.2|l H2.14|l and (|3.1|l we have a real value 



fr'nGwif^hIJif^ 



5h_ 



(3.5) 



The fluctuations around the mean {^.j) are dominated by the noise under our weak signal approximation, and its 
variance erf for the real part of /i^ becomes 



Si{h)Sj{f,) 



Sjj_ 
8Af 



(3.6) 



As the number of Fourier modes Sfi/Af in each segment is much larger than unity, the probability distribution 
function (PDF) for the real part of the measured value /i; is close to Gaussian distribution due to the central limit 
theorem as 



p(Re[/i,]) 



1 



: exp 



2cj 



(3.7) 



Here we neglected the prior information of the spectrum Ogw(/). From eas. ()3.4fl and H3.5|) the signal to noise ratio 
of each segment becomes 



SNRj = 



3Hl 
107r2 



PSi{f)Sj{f) 



Summing up the all the segments quadratically, we get the total signal to noise ratio 



SNR^ 



IOtt^ 



2 \ ^ 



df 



fSi{f)Sj{f) 



(3.8) 



(3.9) 



Note that this expression does not depend on the details of the segmentation {Fi}. We can directly obtain the 
same results by introducing the optimal filter for the product s*j{f)sj{f) to get the highest signal to noise ratio (see 
e.g. H^l). Ea. H3.9|l is given for a single pair {I, J) of the data streams. As the overlap reduction functions 7/j are 
"diagonalized" (7/j ^ only for the following combination {I, J) = (A, A'), {E, E'), (T,T') with / ^ J), the total 
signal to noise ratio given by all of these combinations is evaluated by adding their contributions as 



(SNR) 



2 

opt 



VIOtt^ 



6s 



2E 



df 



fSiif)Sj{f) 



(3.10) 



In this paper we mainly study the case with using all these three combinations, unless otherwise stated. 

Next we discuss how well we can estimate parameters am [fn — 1, Af; M: total number of parameters) that 
characterize the stochastic background spectrum i^icwif)- We can apply standard procedure of the maximum likeli- 
hood analysis for the compressed data {^i} with the probability distribution function (|3.7|l [SI- Then the magnitude 
of the parameter estimation error Aa™ is evaluated by the Fisher information matrix Tmn that is the inverse of the 
error covariance matrix (Ack„iAa„) as 



(AamAa 



/3ff| 



Vl07r2 



obs 



2E 



df 



llj{f?da„MGwif)dc,,nGwif) 

PSiif)Sj{f) 



(3.11) 



The simplest case is the estimation of the amplitude flaw for a flat spectrum ^Icwif) = ^ow = const in the 
frequency range relevant for correlation analysis. In this one parameter estimation with ai ~ IhQgw, the expected 
error becomes 



Aai = 



A17, 



GW 



GW 



{SNRy^. 



(3.12) 



We can easily confirm this by using eas. (|3.10l) and (|3.11() . Actually, the above result (|3.12() holds for the case when 
we estimate only the overall amplitude of the spectrum flaw if) with a known frequency dependence. 
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A more realistic situation is to estimate two parameters, the amplitude ai = In flcw.F and the slope a2 = n, 
assuming that the spectrum has a power-law form 

^Gwif) = nawAf/FT (3.13) 

around a central frequency F. In this case we do not have the simple result H3.12|l for the first parameter ai = In ^lcw,F, 
as the two parameters have correlation. The magnitude of the error Aa2 for the slope does not depend on the choice 
of the frequency F, as we can understand from its geometrical meaning. We can also confirm this directly with using 
eg. (13. 11(1 . In contrast, the error ai and the correlation coefficient r = {Aaia2) /{{Aal') (^Aa^))^^^ depend on the 
frequency F. With a suitable choice of F we can diagonalized the covariance matrix (Aa^aj) {i = 1, 2). Once we get 
the errors Aa; or signal to noise ratio for a specific combination {Uaw,F,n), we can easily obtain the results for a 
different amplitude ^gw,f (but the same slope n) with using scaling relations as 

SNR CX ^GW,F, Atti CX ^GW,F^ ^ ^GW,F- (3-14) 

It is straightforward to confirm these relations with eas. (|3.10|l and (|3.11(l . 

So far we have assumed the situation (weak-signal approximation) that the detector noise {ni{f)} is much larger 
than the gravitational wave background {Hj{f)}. When the latter becomes comparable to the former, the background 
Hj{f) itself becomes a significant source of the fiuctuations for measuring (/i^) in ca. H3.5|l [20l l40l|. For the strong 
signal limit, the signal to noise ratio becomes a asymptotic value of order SNR ^ {TobsfY^'^ that is 10"^ ~ 10^ for 
BBO band with a reasonable observational time. In a recent paper 40] it is shown that for BBO band (> O.lHz) the 
weak signal approximation provides a fairly good estimation for SNR smaller than ^ 200. On the other hand, the 
performance of the Fisher matrix approach becomes worse at low SNR. Therefore, our expressions in this section are 
expected to be reasonable for the background observed with 10 < SNR < 200. 



IV. NUMERICAL RESULTS FOR BBO 



In this section we present numerical results for BBO project. Its primary goal is direct detection of the stochas- 
tic gravitational wave background generated at inflation. As the standard slow-roll inflation predicts a nearly fiat 
spectrum at BBO band, we first assume that the true spectrum has a simple form: ^Gw{f) = ^gw = const. By 
evaluating ea. H3.10|l numerically, we obtain the signal to noise ratio for three possible BBO configurations as 

SNR = im{^^^{^}\^'^ (BBO -lite), (4.1) 

1/2 

SNR = 25.1 ( ) ( ^ ) (BBO - standard), (4.2) 



10-15 J 

f ^GW \ f Fobs 



1/2 

SNR = 251.2 ( ^ ) ( ^ 1 (BBO - grand). (4.3) 



^GW \ I Fobs \ 
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The magnitude ^gw = 10 is close to the upper limit of the amplitude around the BBO band that is consistent 
with the current CMB observation While specification of BBO-lite is not enough to detect the level VLgw ^ 10 
with a sufficient signal to noise ratio, BBO-grand has potential to detect the background close to VIgw ^ 

To calculate the above results we simply integrated ea. (l3.10f) from / = to / = oo. This would be too optimistic 
considering the fact that the frequency below / ^ 0.2Hz might be significantly contaminated by cosmological white 
dwarf binaries p^ . Above ^ 0.2IIz we still have to clean the foreground produced by the binaries made by neutron 
stars or black holes 0, . While it is not clear how well we can actually perform this cleaning (see ^3 for recent 
study), we calculate a less optimistic prediction than eas. H4.1|IH4.2|) and H4.3|l by introducing a lower frequency cut-off 
fcut at fcut = 0.2IIz for the integral H3.10|l . Then the above relations become 



^GW \ ( Toi 



SNR ^ 0.98 ( ^ j (^^J (BBO - lite), (4.4) 

1 /2 

SNR = 17.3 I ( ^] (BBO - standard), (4.5) 



10-16 J \^10yr 

SNR = 116.0 (^) (BBO - grand). (4.6) 
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FIG. 4: Dependence of the signal to noise ratio on the lower cut-off frequency f cut- 
in figure 4 we show how the signal to noise ratio changes with this cut-off frequency. BBO-lite is less sensitive to 
fcut, compared with other two configurations. This is because its minimum of the noise curve is at higher frequency 
than BBO-standard or BBO-grand, as in figure 2. The upper cut-off frequency is not important for our results, if it 
is higher than ~ IHz. Therefore we do not discuss its effects. From figure 2 we can expect that the contribution of 
T-modc to the total SNR is small. Actually, even if we remove (T, T') correlation from eg. 10(1 . the prefactors in 
eas. (|4.4ll(|4.5|) and (|4.3|) change less than 1%. We also study the case with completely aligned two units on a single 
triangle for correlation analysis. The total SNR is obtained by putting 7/j = 56ij and Si{f) = hj{f)^ in ea. H3.10|l . 
Here hi{f) is the effective noise curve for /-mode given in figure 2. We find that the above prefactors become 0.98, 
17.4 and 116.1 respectively, and are very close to the results with the proper transfer functions. Therefore, we can 
approximately discuss performance of a star-like constellation in a very convenient manner with the effective noise 
curves that are often used to represent specification of an interfereometer. 

Now we move to the case with estimating two parameters (ai, = (In flGW,F, n) for the assumed spectrum shape 
^Gwif) = ^GW,F{f / F)" . We fix the true values of the parameters at n = and ^gw,f = 10~^^, and set the lower 
cut-off frequency fcut at 0.2Hz. As we discussed in the last section, the estimation error for the slope n does not 
depend on the choice of the central frequency F ^ and we get 

(Aa^)^^^~ An~ 0.23 (0.039), (4.7) 

for BBO-standard (BBO-grand, respectively). In contrast, the error Aai for the amplitude ^gw,f and the correlation 
coefficient r depend on the frequency F. From observed data we can determine the profile of the spectrum ^Gw{f) 
relatively well around the optimal frequency region where the signal to noise ratio accumulates in eg. H3. 10(1 . However, 
if we take the central frequency F away from this optimal region, the estimated amplitude ^gw,f at the frequency 
F would be strongly affected by the error of the slope n. Consequently, the correlation between errors of the two 
parameters (In ^gw,f^ n) becomes strong. In figure 5 we show 2(T-contour map expected for the two parameter fitting. 
In this figure we take F = IHz that is higher than the optimal frequency region around 0.2 ~ 0.4Hz (see figure 4). 
We can observe a strong correlation between two fitting parameters, and we have \r\ > 0.96 for the results given in 
figure 5. 

We searched the central frequency F that makes the correlation coefficient r = 0, and found F — 0.26Hz (BBO- 
standard) and 0.25Hz (BBO-grand) for the lower cut-off frequency at 0.2Hz. With these choices for the frequency F, 
the error Aai for the amplitude are given by SNR~^ as eg. (|3.12|) . as the two dimensional variance matrix (AaiAaj) 
is now diagonalized. In figure 6 we plot the error ellipses with these frequencies F. Compared with figure 5, the error 
for the amplitude ilGW,F becomes significantly smaller. We can also understand the direction of the ellipse (or the 
sign of the coefficient r) with an argument similar to the discussion given just after ea. (|4.7|l . 

So far we have studied the results mainly for a specific model ^Gwif) = 10~^^ at the BBO band. Here we analyze 
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FIG. 5: Error ellipses {2a; 86%CL) for two dimensional parameter estimation with different lower cut-off frequencies fcut 
(dashed curve: fcut = 0, solid curve: fcut = 0.2Hz, the dotted curve: fcut ~ 0.3Hz). The fisher matrix approach is used. The 
central frequency F is set at _F = IHz that is slightly higher than the optimal sensitivity for measuring the background with 
above two BBO configurations. The true values of the spectrum ^Icwif) are iiGw,F = 10"^*^ and n = 0. 




FIG. 6: Same as figure 5, but with the central frequencies F at 0.26Hz (BBO-standard) and 0.25IIz (BBO-grand). The two 
dimensional covariance matrices are diagonalized with fcut ~ 0.2Hz. 
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FIG. 7: Dependence of various quantities to the slope n. SNR (solid curve) and two errors Aai (dotted curve) and Aa2 
(short-dashed curve) are normalized by their values at n = 0. The long-dashed curve represent the correlation coefficient r. 
The central frequencies are F = 0.26Hz (BBO-standard) and F = 0.25Hz (BBO-grand). 



how various quantities depend on the spectral index n for the two dimensional parameter fitting (ai = IhQgw.f, 
Oi2 = n) with assumed spectral form ^gw = ^GW,F{f /F)'^. We fix ^gw,f — 10~^^, fcut — 0.2Hz, F = 0.26 for 
BBO-standard, 0.25Hz for BBO-grand, but change the slope n. Then the signal to noise ratio, the magnitude of 
errors (Aai, Aa2), and their correlation coefficient r between them are evaluated as functions of the slope n. The 
results are presented in figure 7. Note that the signal to noise ratio and the error for the amplitude ^gw,f depend 
very weakly on the slope. This is because the central frequency F is in the optimal frequency region. With this figure 
and the numerical results given so far, we can get relevant quantities for various combinations {flGW,F,n) by using 
the scaling relations given in ea. H3.14|l . 



V. SUMMARY 



In this paper we studied prospects for the direct measurement of the stochastic gravitational wave background by 
correlation analysis. As a concrete example, we explicitly examined the case with possible BBO configurations that 
use two sets of three spacecrafts to form a star-like constellation. We calculated not only the signal to noise ratio 
for detection, but also how accurately we can measure the basic parameters that characterize the spectrum flGwif): 
such as its amplitude flcw.F or its slope n. 

While it is difficult to detect a level ^Igw ~ 10~^^ with BBO-lite, BBO-grand has potential to detect the background 
close to flGW ~ 10~^^ in 10 years. When we try to measure the amplitude and the slope simultaneously for the spectral 
shape ^Gwif) = ^GW,F{f /Fy\ their errors can be highly correlated and the estimation of the amplitude might be 
degraded, compared with the single parameter fitting for a flat spectrum ^Gw{f) = ^GW- If we take the central 
frequency F around the optimal sensitivity for the correlation analysis, we can nearly diagonalized the two dimensional 
error covariance matrix. 

While we showed the impacts of the low frequency cut-off associated with the potential astrophysical foreground 
that might be difficult to subtract from the data streams, it is not clear how well we can clean the foreground made 
by black hole or neutron star binaries above 0.2IIz. These aspects, especially in relation to the correlation analysis, 
must be clarified to properly understand prospects for the measurement of weak stochastic background {e.g. form 
inflation). In addition we have assumed that the detector noises between combinations {A, A'), {E,E') and {T,T') 
are uncorrelated. In reality they must have correlations to some degree, and estimation of their magnitude is crucial 
for discussing the sensitivity accessible by correlation analysis i2Qj . 
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